Formation and rapid evolution of domain structure at phase transitions in slightly 

inhomogeneous ferroelectrics 



(N 
O 

o 

(N 
C 



o 



A.M. Bratkovsky 1 and A. P. Levanyuk 1 ' 2 
1 Hewlett-Packard Laboratories, 1501 Page Mill Road, Palo Alto, California 94304 
2 Departamento de Fisica de la Materia Condensada, C-III, Universidad Autonoma de Madrid, 28049 Madrid, Spain 

(November 29, 2001) 

We present the first analytical study of stability loss and evolution of domain structure in inhomoge- 
neous ferroelectric samples for exactly solvable model. The model assumes a short-circuited capacitor 
with two regions with slightly different critical temperatures T c i > T C 2, where T c \ — T C 2 <C T c \,Tc%. 
Even a tiny inhomogeneity like 10 _5 K may result in splitting the system into domains below the 
phase transition temperature. At T < T C 2 the domain width a is proportional to (T c i— T)/(T c i— T %) 
and quickly increases with lowering temperature. The minute inhomogeneities in T c may result from 
structural (growth) inhomogeneities which are always present in real samples and a similar role can 
be played by inevitable temperature gradients. 
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The idea that the phase transition in electroded short- 
circuited ferroelectric proceeds into homogeneous mon- 
odomain state [tp is very well known. Similar result 
also applies to free ferroelastic crystals. However, it has 
never been observed. Surprisingly, both electroded fer- 
roelectrics and free ferroelastics do split into domains, 
although they should not. The present paper aims to 
answer why. 

It is generally assumed that in the finite ferroelectric 
samples the domain structure appears in order to reduce 
the depolarizing electric field if there is a nonzero normal 
component of the polarization at the surface of the fer- 
roelectrics [QJ^] (in complete analogy with ferromagnets 
0), if the field cannot be reduced by either conduction 
(usually negligible in ferroelectrics at low temperatures) 
or charge accumulation from environment at the surface 
m. On the other hand, in inhomogeneous ferroelastics 
(e.g. films on a substrate, or inclusions of a new phase 
in a matrix) the elastic domain structure accompanies 
the phase transition in order to minimize the strain en- 
ergy, as is well understood in case of martensitic phase 
transformations || and epitaxial thin films (^-||] . 

In search for reasons of domain appearance in other- 
wise perfect electroded samples, which is not yet under- 
stood, we shall discuss a second order ferroelectric phase 
transition in slightly inhomogeneous electroded sample. 
This problem has not been studied before. We con- 
sider an exactly solvable case of a system, which has two 
slightly different phase transition temperatures in its two 
parts. While the phase transition occurs in the "soft" 
part of the system, the "hard" part may effectively play 
a role of a "dead" layer |l(J and trigger a formation of the 
domain structure in the "soft" part with fringe electric 
fields penetrating the "hard" part. One has to check this 
possibility, but the behavior of the corresponding domain 



structure is expected to be unusual: it should strongly de- 
pend on temperature since further cooling transforms the 
"hard" part into a "soft" one, while the first "soft" part 
becomes "harder" . Since the inhomogeneity is small, one 
might expect that the domains would quickly grow with 
lowering temperature. We indeed find a rapid growth 
of the domain width linearly with temperature in the 
case of slightly inhomogeneous short-circuited ferroelec- 
tric. This behavior is generic and does not depend on 
particular model assumptions. Generally, the inhomo- 
geneous ferroelectric systems pose various fundamental 
problems and currently attract a lot of attention. In par- 
ticular, graded ferroelectric films and ferroelectric super- 
lattices have been shown to have giant pyroelectric pd| ] 
and unusual dielectric response fL2fl , 

We shall consider the case of slightly inhomogeneous 
uniaxial ferroelectric in short-circuited capacitor that 
consists of two layers with slightly different critical tem- 
peratures, so that, for instance, a top part "softens" 
somewhat earlier than the bottom part does. We as- 
sume the easy axis z perpendicular to electrode plates, 
and make use of the Landau free energy functional for 
given potentials on electrodes (zero in the present case) 

§F = F LGD [P] +JdV (-EP - f|) , with 



F LGD [P]= £ [dV[^P* + ^P* 



f(«) 2 + ^Pi], 



(1) 



where P z \ P±\\s the polarization component along (per- 
pendicular to) the "soft" direction, index p — 1(2) marks 
the top (bottom) part of the film: 



< z < h, 



1 



A 2 = A + 5 A, 



-l 2 < z < 0, 



where A%n) — ot(T — T C U2)) and 8 A > (meaning T c2 < 
T c i). The constant a = 1/Tq, where To ~ T at (~ T c ) for 
displacive (order-disorder) type ferroelectrics. 

The equation of state is SF LGD [P]/SP = E = — Vy>, 
where tp is the electrostatic potential, or in both parts of 
the filmp= 1,2 : 

E z = -O z if = A p P z + BP* - DV\P Z - gd 2 P z , (2) 
E X =A 1 _P X , (3) 

These equations should be solved together with the 
Maxwell equation, div(E + 4ttP) = 0, or 



(d 2 z + e a Wl) <p = 4nd z P z 



(4) 



where the dielectric constant in the plane of the film is 
e a = l + 4ir/A ± . 

Loss of stability.— We shall now find conditions for loss 
of stability of the paraelectric phase close to T c \ with 
respect to inhomogeneous polarization. At the point of 
stability loss the polarization is small and the nonlinear 
term P z must be omitted. We are looking for a nontrivial 
solution in a form of the " polarization wave" , 



P z ,ip(x e 1 ' 



(5) 



We shall check later that the stability will be lost for the 
wave vector kl\ > 1, so that V\P Z = k 2 P z > gd 2 P z ~ 
P z /l 2 , and the last term in the right-hand side of (§) 
should be dropped. Going over to Fourier harmonics in- 
dicated by the subscript k, we obtain 



9k 



e a k 2 (fi k = ^P' zk , 



(6) 



where the prime indicates derivative ( /' = df/dz, f" = 
d 2 f/dz 2 ). We can exclude P z k with the use of the lin- 
earized equation of state (0), which gives 



(A p + Dk 2 )P; 



zk • 



(7) 



fl - £ L Vk = 0, where 



Substituting 

this into (|^), we obtain if" 

we have used \A + Dk 2 \ /4tt <C 1, which is always valid 
in ferroelectrics. We shall see momentarily that the non- 
trivial solution appears only when A\ + Dk 2 < 0, while 
A2 + Dk 2 > 0. The resulting system is 



<Pik + Xi* Vk = 0, 
<P2k ~ xlk 2 ^2k = 0, 



where x? = j r 

ary condition reads as 



s xi 



- t (A 2 +Dk 2 
4tt 



<P2k 



A 1 + Dk 2 A 2 + Dk 2 ' 



(8) 
(9) 

The bound- 



(10) 



where we have used \Ai + Dk 2 1 /An <C 1. We obtain from 
Eqs. (pj|)-(|Tc|) the condition for a nontrivial solution 



XitanxifcZi = X2tanhx2fc^2, 



(11) 



which has a homogeneous solution k = and the in- 
homogeneous solution with k — k c (|l3|), hence we have 
to determine which one is actually realized. The inho- 
mogeneous solution is easily found for X2M2 1, where 
tanh can be replaced by unity. Close to the transition 
X2/X1 ^ 1) and the solution is 



Xikh = 



7T X2M1 



2 1+xafcJi 



7T 

2' 



x k 2 i- 



(12) 



.There is 



when X2k.l1 > 1. This gives \A\ = Dk 2 

no solution for \\ < 0. The minimal value of A for the 
nontrivial solution (onset of instability) is defined by 



e a Dl 2 
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\A\ C = 2Dk 2 = 



27r V2 D l/2 2 ^3/2 , 



1/2, 



M 2 h ' 



(13) 
(14) 



where we have introduced the "atomic" size d a t ~ vD 
comparable to the lattice parameter. We obtain the cor- 
responding tiny shift in the critical temperature [see es- 
timates below Eq.(|l6|)] T& — T c ~ Tod a t/eV 2 h- Hence, 
the system looses its stability very quickly below the 
bulk transition temperature. It is readily checked that 
the assumptions we used to obtain the solution are eas- 
ily satisfied. Indeed, X2^2 ^ 1 and X2kli ^S> 1 both 
correspond to approximately the same condition when 
l\ ~ I2 '■ SA ^> — 4 1/2 meaning that the difference 

between T c should be larger than the shift of T c . 

Now we have to determine when the transition into 
inhomogeneous state occurs prior to a loss of stability 
with respect to a homogeneous polarization. The homo- 
geneous loss of stability corresponds to A = Ah found 
from 



A h h + {A h + SA) h = 0. 



(15) 



For the inhomogeneous state to appear first, there must 

be A c > A h , or SA > ^Mp^- ^. This means that 

very tiny inhomogeneity in the sample is enough to split 
it into the domain structure, 



T c i — T C 2 — Tq 



^/ 2 {h+l 2 ) d at 



(16) 



which is estimated as T at < e a (10 4 — 

lO^lCT 7 ^ = (10^ 3 - 10" 2 )K Q for displacive systems, 
and T c 4fi T - < (l0~ 5 - 10" 4 )K for order-disorder sys- 
tems. Certainly, such a small temperature and/or com- 
positional inhomogeneity exists in all usual experiments. 
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with the boundary conditions at the interface z = 



(22) 



The corresponding electrostatic (stray) field part of the 
energy is found as |l0| F es = | J dAcr s (p (z = 0) , where 
<7 S is the density of bound charge at the interface , corre- 
sponding to only the spontaneous part of the polarization 
Pqi (x), and integration goes over the area A between two 
parts of the film. We calculate this expression by going 
over to Fourier expansion ( |l9| ) and using the fact that in 
the present geometry a s (x) = — Poi(x) (and, therefore, 
its Fourier component a s k = —Poik), 



FIG. 1. Schematic of the domain structure with the pe- 
riod 2a in inhomogeneous ferroelectric film of the thickness 
h + l 2 . Top and bottom layers have slightly different critical 
temperatures T c ± > T c2 , T cl - T c2 < T cl ,T c2 . Slightly be- 
low T c i the top layer splits into domains with electric fringe 
field propagating into the bottom layer (fringe field shown as 
the hatched area in the top panel). The domains persist and 
evolve below T c2 when both layers exhibit a ferroelectric (or 
ferroelastic) transition (bottom panel). 



Domain structure at T C 2 < T < T c \ (A < 0, A + 5 A > 
0).— After stability loss the resulting "polarization wave" 
quickly develops into a domain structure, as we shall now 
demonstrate. In the region well below T c i we can use the 
linearized equation of state 

E z = (A + SBP^)(P Z - Pqi) = -2A(P Z - P i), (17) 



where |pn| = \/—A/B is the spontaneous polarization 
in the top layer, which gives P zl = P i + wm^zi P Z 2 = 
■j^E z , for the top and bottom layers, respectively. In 
this case the equation for the potential ip (0) reduces to 
a standard Laplace equation (e c d z + e a V j_J f = 0, with 
the boundary condition 



e c id z (pi - e c2 d z ip 2 = 47rP i(x), 



(18) 



where e cl = 1 + 2n/\A\, e c2 = 1 + -i7r/A 2 . 

The spontaneous polarization in the top layer alter- 
nates from domain to domain as Poi(x) = ±|Poi| = 
±y/—A/B. We are looking for a solution in a form of 
a domain structure with a period T — 2a (Fig. |l|), 



Poi(x) =Y, P oike lkx , <p[x) = 



k 



(19) 



with k = k n — InnjT — irn/a, n = ±1,±2,... Going 
over to the Fourier harmonics, we can write the Laplace 
equations for both parts of the film as 



£ak 2 <fi2k = 0, 



£ c2^2fc 



(20) 
(21) 
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~A 



y-' 47r|P i fc | 

fe>0 



kD k 



(23) 



n 1/2 



el{ 2 coth .l—kli + eto* coth <l—kl 2 



1/2, 



, with 

k = nn/a, n = 1,2,..., similar to jL3j. Note that here 
Poifc = 2 |Poi| /inn, n — 2j + 1, j = 0, 1, ... and zero oth- 
erwise. Adding the surface energy of the domain walls, 
we obtain the free energy of the domain pattern 
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(24) 



where D n = Dp. . Not very close to T c the argument of 



■c2 



coth is j^kli > 1, so that Dk 
Minimizing the free energy, we find the domain width 
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14C(3) 



AiZi 



1/2 



(25) 



where Ai = Ti/Pqi = datl^l 1 ^ 2 is the characteristic mi- 
croscopic length, and d at = 2-^-D 1 / 2 is comparable to 
a lattice spacing ( "atomic" length scale) . The expression 

(H) is valid when Jj^kh > 1, or \A\ > 2d at / (wel^l^ , 
meaning that one has to be below T c by a tiny amount 
T c i — T > Tod at / (el/ 2 li) , estimated earlier. Note that 
close to T c i one obtains for the domain width 



a K 



^ 2 el /2 
7\/2C(3) 



d a th 



1/2 



(26) 



and this value does not depend on temperature. We shall 
formally refer to this result as the Kittel domain width. 
Incidentally, close to T c2 the domain width is a « 

1/2 



tM /2 4 / 2 

new 



-1/4 



e c2 4 ; which formally diverges oc 



(T — T c2 ) 1/-± . However, in the vicinity of T c2 the in- 
duced polarization in the formerly "hard" part has about 
the same value as the spontaneous polarization in the 
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"soft" part, P Z 2 ~ Pqi- Then the equation of state 
in the bottom part becomes strongly non-linear, since 
the cubic term is much larger than the linear term, 
BPl 2 w BPl x = AP 01 w AP z2 > A 2 P z2l in the equation 
of state (since A ^> A 2 close to T c2 ), so the response of 
the bottom layer does not actually soften in this region. 
In this case our derivation does not apply, but it is prac- 
tically certain that the domain structure in the vicinity 
of T c2 would evolve continuously upon cooling, Fig. 0. 

Domain structure at low temperatures (T < T c2 , 
A < 0, A + 5 A < 0).- When the system is cooled to 
below the critical temperature T c2l a spontaneous po- 
larization IP02I = V / —A 2 /B also appears in the bottom 
layer. The domain structure simultaneously develops in 
the whole crystal with domain walls running parallel to 
the ferroelectric axis through the whole crystal (if they 
were discontinuous at the interface between the two parts 
of the crystal this would have created a large depolarizing 
electric field). The electrostatic energy requires a solu- 
tion of the same Laplace equations ( pp[ ) and (^l]), only 
the boundary condition (E2h would now read 



(27) 



e cl<fi'lk - £c2<P 2 k = 47r(P ifc - Po2k), 



where e c i( 2 ) = 1 + 27r/|Ai(2)| ~ 27r/| J 4 1 ( 2 ) I- Note that the 
density of the bound charge at the interface, correspond- 
ing to this discontinuity of spontaneous polarization, is 
now Ok = — (Poik — Po2k)- Therefore, we immediately ob- 
tain for the total free energy of the structure, analogously 
to the previous case (pjf). 



F = P&Axh + Pj 2 &2h 
A a 

, 16(Pqi - P 02 ) 2 a 



E 



" {2 3 + iyD 2j+1 



(28) 
(29) 



where Am 



(2) 



dat^J\A 



1(2)1- 



Not very close to T c2 we 
would have </—kl 2 > 1, and the minimum of the free 
energy F is achieved for the domain width 
1 



1 



- P02/P01 
2 1/2 / 1/2 



£ l/2 
e c2 



14C(3) 



A1Z1 



p2 



1/2 



(30) 



Close to the critical point T c2 the domain width formally 
behaves as a oc e^ 4 oc (T c2 — T) -1 / 4 , as found just above 
T c2 before. The same argument indicates though that 
our derivation does not apply in this region, but non- 
linearity should not cause a substantial change in the 
domain structure. 

With lowering the temperature to the region where 
\A\ > 5 A, we will have P02/P01 = yJ{A + 5 A) /A w 1 + 
5A/2A, so that 1 - P02/P01 ~ 2\A\/5A > 1 becomes 



a/a 



cz cl 



FIG. 2. The domain width in slightly inhomogeneous fer- 
roelectric or ferroelastic in the units of an, the Kittel width 
(^6j). a — ax when the domain structure sets in at T w T c i, 
and then it grows linearly with the temperature to large val- 
ues a 3> aK ■ 



a large prefactor. Note that in this region e c ± w e c i = 
2tt/\A\, Ai w A 2 = d a ty/\A\, and the domain width 
evolves as 



6A 



2 5/2^5/2 e l/2 

7C(3) 



d at (h + h) 



1/2 



(31) 



It becomes much larger than the Kittel width, 

» 1, (32) 



= 2 3/2 f' h + h \ T-\ T 



OA" 



<:2 



growing linearly with lowering temperature, if the pin- 
ning of the domain walls is neglig ible (Fig. |). Close to 
the lower critical point the linearized equation of state 
does not apply but the response of the bottom layer re- 
mains finite, and we expect, as mentioned above, that the 
domain structure would evolve rather gradually across 
T C 2, Fig. I 

Summarizing, in a ferroelectric sample with a tiny in- 
homogeneity of either the critical temperature or tem- 
perature itself (i.e. in the presence of a slight tempera- 
ture gradient and/or minute compositional inhomogene- 
ity across the system) the domain structure abruptly 
sets in when the spontaneous polarization appears in the 
softest part of the sample (i.e. the part with maximal 
T c ). This takes place when the difference in T c in the 
parts of the sample is just (1CP 3 — 10 _2 )K for displacive 
systems, and even smaller, (lO -5 — 10~ 4 )K, for order- 
disorder systems. The period of the structure then grows 
linearly with lowering temperature and quickly becomes 
much larger than the corresponding Kittel period. 

This result does not depend on specific geometry as- 
sumed in the present model. Indeed, if local T c — T c (z) 
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varies continuously, like in graded ferroelectrics [Ly], it 
can be approximated by a piece-wise distribution of a se- 
quence of "slices" . Upon cooling the system first looses 
stability in the softest part of thickness l s , which is de- 
rived from the position of the boundary where local 
T c = 0, with respect to a domain structure with fine 
period oc \fl~ s . The domains extend into the bulk of the 
system and become wider with further cooling, since l s 
increases. In electroded sample there will be no branch- 
ing and domain walls would run straight across all trans- 
formed slices. Otherwise, discontinuities would have re- 
sulted in very strong depolarizing field. If the overall 
inhomogeneity is small, the picture would obviously re- 
main very similar to the two-slice model solved above. 
The same arguments remain valid if the inhomogeneity 
were to have more complex form/distribution in a sam- 
ple. The novel feature of the present effect of the depolar- 
izing field is that it appears not due to surface charges, 
which are screened out by the electrodes, but because 
of the bulk inhomogeneity. The bulk depolarizing fields 
are present in other important classes of inhomogeneous 
ferroelectrics, graded ferroelectric films ]iT| and super- 
lattices of different ferroelectrics (e.g. KNbOa/KTaOa) 
]l2fl , and may be responsible for their unusual behavior. 

We have shown that a very tiny temperature gradi- 
ent, or a slight compositional inhomogeneity, etc., would 
result in practically any crystal eventually splitting into 
domains no matter how high the quality of it is. The 
unusual evolution of the domain pattern, found in the 
present paper, when it starts from very fine domains 
at T c and then grows linearly with temperature to very 
large sizes, has been reported in Ref. |L4| for ~ 1mm 
thick TGS crystals. It is worth noting that the result is 
very general and applies also to slightly inhomogeneous 
free ferroelastic crystals Ju|. Other implications include 
extensively studied graded films and ferroelectric super- 
lattices [[y]|lj]. It would be very interesting to perform 
controlled experiments for the domain structure close to 
the phase transition to check the present theory. 
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